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Abstract 

In this talk we discuss the analogy between data from heavy-ion collisions and the Cosmic Microwave Back- 
ground. We identify p T correlations data as the heavy-ion analogy to the CMB and extract a power-spectrum from the 
heavy-ion data. We define the ratio of the final state power-spectrum to the initial coordinate-space eccentricity as the 
transfer-function. From the transfer-function we find that higher n terms are suppressed and we argue that the suppres- 
sion provides information on length scales like the mean-free-path. We make a rough estimate of the mean-free-path 
and find that it is larger than estimates based on the centrality dependence of V2. 
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A commonly quoted goal of the heavy-ion programs at Brookhaven National Laboratory and CERN is to recreate 
conditions similar to those shortly after the Big Bang [ \\. Learning about the early stages of the matter produced in 
heavy-ion collisions from the observation of hadrons in the final state is in some sense similar to understanding the 
stages of the early universe from the observation of the Cosmic Microwave Background (CMB) |2|. In this talk we 
explore the analogy between heavy-ion collisions and Big Bang cosmology 0. In particular, we present the heavy- 
ion equivalent of the CMB and from that we determine the power-spectrum for the little bangs. We then estimate the 
transfer-function necessary to produce the spectrum from the initial conditions of the collisions. 

Quantum fluctuations from the early universe show up as hotspots at the surface of last scattering, creating the 
CMB measured for example by WMAP. Measurements of the CMB reveal temperature fluctuations corresponding 
to over- or under-densities present at the surface of last scattering at about 400,000 years after the Big Bang |2|. 
These density fluctuations ultimately explain the structure in our universe. Maps of the temperature of the CMB are 
determined from measurements of the blackbody spectra at 2M points in the sky. The temperature is found to be 
very smooth with relative variations only showing up at approximately 10~ 5 . The fluctuations tend to only be at short 
distances; small lumps not large ones. This result has been explained by an inflationary period when the larger scale 
fluctuations were pushed outside of the horizon, causally separating them. 

The scale of the correlations can be most easily studied by extracting a power-spectrum from the CMB. The 
power-spectrum shows that most power is for large wave-numbers (meaning small wavelength). The lack of power at 
small wave-number is strong evidence of inflation. Peaks in the power-spectrum are caused by acoustic phenomena in 
the early universe as density perturbations in the universe propagate as sound waves. This gives rise to the structure in 
the anisotropies of the microwave background, notably a characteristic angular scale and the famous acoustic peaks. 
From this spectrum the acoustic peaks are fit with models to extract cosmological parameters. In this talk we explore 
an analogous measurement in heavy-ion collisions. We determine the power-spectrum for the little bangs and we 
estimate the transfer-function necessary to produce the spectrum from the initial conditions. 

Just as quantum fluctuations stretched to cosmic sizes by inflation show up in the CMB, we expect fluctuations 
from the beginning of the little bangs to show up in heavy-ion data. Mishra et. al first have explored the analogy 
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Figure 1: Top: the length scales in the CMB. Bottom: the length scales probed with higher n in heavy-ion collisions. 

between the expansion of heavy-ion collisions starting from a lumpy initial energy density and the expansion of the 
universe starting with quantum fluctuations stretched to cosmological sizes, and based on that analogy suggested mea- 
suring RMS values of the Fourier coefficients v„ (known as flow coefficients) as the power-spectrum Q. They didn't, 
however, make a connection between the RMS of v„ and the already existing two-particle correlation measurements. 
Sorensen later pointed out [4 1 that many of the surprising correlation structures seen in RHIC data (the ridge and mach 
cones [5 1) may be understood in terms of fluctuations of higher Fourier components of v„ which Mishra et. al. had 
argued would arise from anisotropics in the initial energy density converted into momentum-space during the expan- 
sion. Several other authors have also previously explored the connection between initial conditions and two-particle 
correlations @ M CED - 

While quantum fluctuations stretched to cosmic sizes by inflation explain the temperature fluctuations in the 
CMB, lets see how density fluctuations at the beginning of heavy-ion collisions might show up in the correlations 
data. We have some expectations about the initial state: given that the collisions involve nuclei with a finite number 
of nucleons, we expect that the initial density is inhomogeneous, it has a lumpiness. Every event is unique and no 
event is homogeneous. We do not dwell for the moment on the exact theory to describe these initial fluctuations, 
but only point out that the finite number of charges in the approaching nuclei implies that the initial (energy) density 
is lumpy (even if we consider a larger number of gluons in the nucleus, those gluons are still confined in the radial 
direction within the typical size of the nucleons). The lumpy initial density may survive through the evolution and 
show up at the end as hotspots in the surface of last scattering. We expect this because the rapid expansion boosts 
these lumps outward and the decay of the hotspots in coordinate-space manifests in correlations in momentum-space 
due to flow [6]. The lumpy initial conditions show up in the 2-particle correlation structures. The plausibility of this 
conjecture has been verified with multiple models for heavy-ion collisions ||9l [TOl ITTI [T2l . 

For the analogy with the CMB we are interested in temperature-correlations. Thus we need to consider how 
much of the initial inhomogeneity will be transferred into the final state. Simulations based on hydrodynamic models 
indicate that collisions of Au nuclei (12 fm across) may contain hotspots of size 1.5 fm and that remnants of those 
hotspots persists during the collisions evolution [12J. How much of this initial inhomogeneity is being transferred 
into the final state depends on the mean-free-path l m f p and it's relation to the size of the hotspots. The higher the 
value of n, the smaller the length scale we are probing. An n — 2 asymmetry (think of V2) can be translated because 
the l m f p is small compared to the length probed by n = 2. That's why ideal hydro works well for Vz- The fact that 
hydrodynamic models do a reasonable job of predicting the value of vo suggests that l m f p can be considered small 
compared to the size of the system. But what about n — 3 or even larger? As illustrated in Fig. [T] as we increase n, 
we reduce the length scale probed. If the mean-free-path is greater than the length scale, there is no translation from 
coordinate-space into momentum-space. We only expect an efficient conversion of coordinate-space anisotropies into 
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momentum-space when l m f p < 2n(R}/n where (R) is the average radial position of the systems constituents. This 
relation was checked in a simulation by calculating the initial participant eccentricity for all values of n part ) ifTUl 
from a Monte Carlo Glauber model lfl3l . For central collisions, when the participants are treated as point like, e\ par , 
is nearly independent of n. When the participants are smeared over a region of size r part , e\ pan at higher n becomes 
quenched. We find that e\ part is reduced by half when r pari < 2n(R)/n and argue that l m f p will have a similar effect as 
r pan- The power-spectrum tells us the power transferred into each harmonic n as a function of n. Clearly, the behavior 
of the power-spectrum can be used to extract information about length scales like l m / p . By examining the power- 
spectrum of heavy-ion collisions which includes information for all values of n (beyond just n — 2 or 3), we hope to 
better constrain l m j p [14|. We note that other length scales that could be important have been identified: thermal scale; 
time of free streaming; size of hadrons; etc 1 15 1. 

To determine the power-spectrum we take the STAR data of p T correlations vs the relative azimuthal angles and 
the relative rapidity between the particles lfT6l . We use this correlation measurement since this is a good proxy for 
point-to-point temperature-fluctuations: just like the blackbody spectrum, the pj spectra can reflect the temperature in 
heavy-ion collisions. A narrow peak positioned around small angle separation is observed in the data. This tells us that 
if a particle comes out with above average pj, then the nearby particles also tend to have large pj . This is consistent 
with expectations from hotspots on the surface of last scattering. The correlation of these fast particles suggests that 
they are born out of the same high-density, high-temperature lump, i.e. the same hotspot. The relevance of temperature 
fluctuations to heavy-ion data is supported by the improvement to pj spectra fits using a Tsallis distribution [ 17 1. 

We considered a slice of the pj correlations at zero relative rapidity vs relative azimuthal angle. We Fourier- 
transform this correlation to generate a power-spectrum vs harmonic number n. The coefficients of the Fourier de- 
composition a n give us the power-spectrum. This is shown in the left plot of Figure [2] The a n are similar to v\ but 
extracted from p T correlations instead of number correlations IfTSl . a n is the power-spectrum for heavy-ion collisions 
and it is the analogous measurement to that of the CMB. We find that most of the power is in the lowest modes and a„ 
falls off rapidly with n. Modes other than n = 2 however still show a significant contribution. 

For the CMB, the lowest modes are strongly suppressed and significant acoustic oscillations show up at I — 200. 
This is due to the fact that in the early universe inflation stretches the fluctuations by 30 orders of magnitude so that 
the small wave numbers are outside the acoustic horizon. Mishra et. al had proposed looking for suppression of small 
n modes due to an acoustic horizon in heavy-ion collisions 0, but we see is clearly a suppression of large n modes. 
We argue that this is because viscosity is far more important in heavy-ion collisions than in the early universe lfT9ll . 

Apparently, the damping due to various length-scales in heavy-ion collisions (l m f p , etc.) is the dominant effect. 
Next we want to compare the final power-spectrum to the initial coordinate-space anisotropy. The azimuthal distribu- 
tion of matter in the initial overlap zone is calculated using the participant eccentricity for all harmonics n s\ parl [ 10 1. 
We start with a Glauber model 1 13 1 for the initial spatial distribution of participants for perfectly central collision (b=0 
fm) in order to compare to data also for the 5 % most central collisions. The results are shown in the middle plot of 
Figure [2] The large n = 2 term persists even for central collisions because we include the intrinsic deformation of the 
Au nucleus in our Monte Carlo. The n — 1 term is small because the participants are re-centered to the center of mass. 
For central collisions for n 1 2 the eccentricity is nearly independent of n (the higher harmonics are important too). 
This is because for symmetric collisions and point-like participants, e\ parr depends only on the number of participants, 
independent of n. We expect, however, that the conversion of higher harmonic eccentricity will be damped due to the 
existence of the length scale l m f p . Giving the participants a size causes the eccentricity to fall off faster with n. The 
drop occurs where the size of the participant is larger than the lengthscale corresponding to n the s\ part is smeared 

out (20). 

We define the ratio a n /sl part as the transfer-function [14] and show it on the right plot of Figure^ The transfer- 
function gives us information about how efficiently the coordinate-space anisotropy has been convertedinto momentum- 
space. The transfer-function derived from data and the Glauber model shows that as n increases, the efficiency for 
converting coordinate-space anisotropy into momentum-space quickly drops off. This can be expected from the 
condition that the transfer-function should go to zero when l m f p is larger than the length-scale associated with n: 
n > 2n{R)ll m r p . We can make a crude estimate for l m f p based on the transfer-function. If we take (R) = 3 fm for the 
average radial position of a participant and n = 6 as the harmonic beyond which the conversion is inefficient, then we 
get l m f p = 3 fm. This estimate can be improved with a more complete model of heavy-ion collisions. This estimate 
corresponds to a viscosity several times larger than estimates based on the centrality dependence of V2 lfl9l or from 
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Figure 2: Left: the power-spectrum a„ extracted from pj correlations in heavy-ion collisions. Middle: the participant eccentricity calculated from 
a Monte-Carlo Glauber model. Right: the transfer-function for heavy-ion collisions showing the suppression at higher values of n in the final 
momentum-space compared to the initial coordinate-space. 

the longitudinal width of pj correlations ETI . Our crude estimate, however, is geometry based, not accounting for 
the various phases of the expansion. 

We've used p T correlations from heavy-ion collisions to extract a power-spectrum analogous to the power- 
spectrum extracted from the Cosmic Microwave Background. We calculated the initial state eccentricity and then 
find the transfer-function which we define as the ratio of the power-spectrum to the initial eccentricity. We argue 
that the transfer-function required to describe the RHIC data may potentially be understood in terms of an ineffi- 
ciency in conversion of coordinate-space anisotropy into momentum-space when l m f p > 2n{R)ln. We used this 
transfer-function to make a rough estimate of the mean-free-path of the system's constituents. The introduction of 
the harmonic axis n provides more detailed information. This approach represents a new method for determining the 
characteristics of heavy-ion collisions and the QGP and should be further investigated. It stands as a challenge for 
models of heavy-ion collisions to correctly describe this transfer-function: in other words, to correctly describe the 
conversion of coordinate-space anisotropies at all scales into momentum-space anisotropies at all scales. 
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